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Abstract. In 1934, Jordan et al. gave a necessary algebraic condition, the Jordan identity, for 
a sensible theory of quantum mechanics. All but one of the algebras that satisfy this condition 
can be described by Hermitian matrices over the complexes or quaternions. The remaining, 
exceptional Jordan algebra can be described by 3 x 3 Hermitian matrices over the octonions. 

We first review properties of the octonions and the exceptional Jordan algebra, including 
our previous work on the octonionic Jordan eigenvalue problem. We then examine a particular 
real, noncompact form of the Lie group E^, which preserves determinants in the exceptional 
Jordan algebra. 

Finally, we describe a possible symmetry-breaking scenario within E§: first choose one of 
the octonionic directions to be special, then choose one of the 2x2 submatrices inside the 3x3 
matrices to be special. Making only these two choices, we are able to describe many properties 
of leptons in a natural way. We further speculate on the ways in which quarks might be similarly 
encoded. 



1. Introduction 

A personal note from Corinne: During the academic year 1986/87, Tevian and I were living in 
York, newly married and young postdocs. Tevian was working in the mathematics department 
there, doing research in general relativity, and I was working in Durham, with David Fairlie, 
just beginning my research into the octonionic structures reported here. Imagine my pleasure, 
when I found out that York had its own resident expert on the octonions! I returned to York 
the following summer, to work with Tony on an attempt to describe the superstring using 
octonions [1]. I will be forever grateful to him, not only for the generous way in which he shared 
his vast knowledge and experience in this field, but also for the friendship, respect, collegiality, 
and mentorship, which he also generously shared. 



2. Exceptional Quantum Mechanics 

In the Dirac formulation of quantum mechanics, a quantum mechanical state is represented by 
a complex vector v, often written as \v), which is usually normalized such that v'v = 1. In the 
Jordan formulation [2, 3, 4, 5], the same state is instead represented by the Hermitian matrix 
vv\ also written as which squares to itself and has trace 1. The matrix vv^ is thus the 

projection operator for the state v, which can also be viewed as a pure state in the density 
matrix formulation of quantum mechanics. Note that the usual phase freedom in v is no longer 
present in vv\ which is uniquely determined by the state (and the normalization condition). 

A fundamental object in the Dirac formalism is the probability amplitude v^w, or (v\w), 
which is not however measurable; it is the squared norm |(t;|w;)| 2 = (v\w)(w\v) of the probability 



amplitude which yields measurable probabilities. One of the basic observations which leads to 
the Jordan formalism is that these probabilities can be expressed entirely in terms of the Jordan 
product of projection operators, since 

(v\w){w\v) = (v^w)(w^v) = tr(vv^ o ww^) (1) 

where o denotes the Jordan product [2, 3] 

A°B = ^(AB + BA) (2) 

which is commutative but not associative. 

Remarkably, the Jordan formulation of quantum mechanics does not require (v and) A to be 
complex, but only that the Jordan identity 

(AoB)oA 2 = Ao(BoA 2 ) (3) 

hold for two Hermitian matrices A and B. As shown in [3], the Jordan identity (3) is equivalent 
to power associativity, which ensures that arbitrary powers of Jordan matrices — and hence of 
quantum mechanical observables — are well-defined. 

The Jordan identity (3) is the defining property of a Jordan algebra [2], and is clearly satisfied 
if the operator algebra is associative, which will be the case if the elements of the Hermitian 
matrices A, B themselves lie in an associative algebra. Remarkably, one further possibility 
exists, for which the elements of the Hermitian matrices do not lie in an associative algebra. 
This example is the Albert algebra (also called the exceptional Jordan algebra) H^(Q) of 3 x 3 
octonionic Hermitian matrices [3, 6]. 1 In what follows we will restrict our attention to this 
exceptional case. 2 



3. Quaternions and Octonions 

The Hurwitz Theorem states that the real numbers M, complexes C, quaternions H, and 
octonions O are the only (normed) division algebras (over the real numbers). 3 The quaternions 
and octonions are extensions of the familiar real and complex numbers. A quaternion is an 
arbitrary real linear combination of the real identity element 1 and three different square 
roots of minus one, which are conventionally called {i,j, k} and satisfy the multiplication table 
given in Figure 1. Similarly, the octonions are formed from seven square roots of minus one 
which we will call {i,j,k,k£,j£,il,£}, whose multiplication table is summarized in Figure 2. In 
these multiplication tables, each point corresponds to an imaginary unit. Each line or circle 
corresponds to a quaternionic triple with the arrow giving the orientation. For example, 

kl = kl (4) 
iki = k (5) 
klk = I (6) 

and each of these products anticommutes, that is, reversing the order contributes a minus sign. 

1 The 2x2 octonionic Hermitian matrices ^(O) also form a Jordan algebra, but, even though the octonions 
are not associative, it is possible to find an associative algebra that leads to the same Jordan algebra [3, 7]. 

2 In this case, the equivalence in (1) fails; it is the right-hand side which provides the correct generalization. 

3 A division algebra is a vector space over a field (in this case R) which is also a ring with identity under 
multiplication, and in which ax — b can be uniquely solved for x (unless a — 0). A normed division algebra 
satisfies (12) in addition, and is therefore also an integral domain, that is, a ring in which ab — implies a = or 
b = 0. 
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Figure 1. The quaternionic multiplication table. 
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Figure 2. The octonionic multiplication table. Each line segment should be thought of as 
circle, identical to the quaternionic multiplication table in Figure 1. 

We define the conjugate a of a quaternion or octonion a as the (real) linear map which reverses 
the sign of each imaginary unit. Thus, if 

a = a\ + d2 i + 03 j + ct4 k + 05 k£ + a§ j£ + a-j il + as i (7) 

then 

a = cti — (Z2 i — 03 j ' — a4 k — 05 k£ — a& j£ — a-j ii — a% I (8) 
Direct computation shows that 

ab = ba (9) 
The norm \a\ of an octonion a is defined by 

|o| 2 = aa = af + a\ + a\ + a\ + a 5 + a\ + ay + a\ (10) 



The only octonion with norm is 0, and every nonzero octonion has a unique inverse, namely 



a" 1 = -p-^ (11) 



For all the normed division algebras, the norm satisfies the identity 

|o6| = |a||6| (12) 

A remarkable property of the octonions is that they are not associative! For example, compare 

(ij)(£) = +(k)(t) = +ki (13) 
(i)(j£) = (i)(j£) = -k£ (14) 

However, the octonions are alternative, that is, products involving no more than 2 independent 
octonions do associate. The commutator of two octonions a, b is given as usual by 

[a, b] = ab — ba (15) 

and we define the associator of three octonions a, b, c by 

[a, b, c] = (ab)c — a(bc) (16) 

which quantifies the lack of associativity. More generally, both the commutator and associator 
are antisymmetric, that is, interchanging any two arguments changes the result by a minus 
sign; replacing any argument by its conjugate has the same effect, because the real parts don't 
contribute to the associator. 

The units i, j, k, k£, j£, i£, and £ are by no means the only square roots of —1. Rather, 
any pure imaginary quaternion or octonion squares to a negative number, so it is only necessary 
to choose its norm to be 1 in order to get a square root of —1. The imaginary quaternions of 
norm 1 form a 2-sphere in the 3-dimensional space of imaginary quaternions. The imaginary 
octonions of norm 1 form a 6-sphere in the 7-dimensional space of imaginary octonions. 

Any such unit imaginary quaternion or octonion s can be used to construct a complex 
subalgebra of H or O, which we will also denote by C, and which takes the form 

C = {a R + a s s} (17) 

with or, a s £ ML Regarding s as the complex unit, we have the familiar Euler identity 

e se = cos 9 + s sin 9 (18) 

so that any quaternion or octonion can be written in the form 

a = re w (19) 

where 

r = \a\ (20) 

Any two unit imaginary octonions s and t that point in independent directions determine a 
quatemionic subalgebra of O. 
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Figure 3. One class of elements of G2. These transformations are also contained in the preferred 
517(3) that fixes t. 

4. The Structure of G 2 and SU(3) 

The freedom to choose an entire 2-sphere or 6-sphere of square roots of minus one within 
the 3-dimensional space of the pure imaginary quaternions or the 7-dimensional space of pure 
imaginary octonions leads one to investigate the transformations that preserve the corresponding 
multiplication table. These transformations form the automorphism group of the corresponding 
division algebra. 

In the case of the quaternions, one can imagine rotating i to any pure imaginary point on the 
2-sphere (2 degrees of freedom). Then j can be chosen to be any direction perpendicular to the 
direction of i, i.e. on the equator of the resulting 2-sphere (1 degree of freedom). The direction 
of k is determined by the multiplication table. The 3-dimensional automorphism group of the 
quaternions is therefore seen to be 50(3). 

For the octonions, one can again imagine rotating i to any pure imaginary point on the 
6-sphere (6 degrees of freedom). Then j must again be perpendicular to i (5 degrees of freedom) 
and the direction of k is fixed by the multiplication table. But I is now free to be any direction 
perpendicular to all of the i, j, and k directions (3 degrees of freedom) and the directions of 
the remaining units are determined by the multiplication table. This 14-dimensional Lie group 
turns out to be the exceptional group G2. 

Another way of envisioning the transformations in the group G2 was first shown to us by 
Sudbery [8]. Consider the octonionic unit kt at the top of the multiplication table shown in 
Figure 3. There are three pairs of octonionic units that form quaternionic subalgebras with 
kt, i.e. {j,it}, {jt, i}, and {k,£}. We call these the pairs that "point to" kt. If the elements 
of two of these pairs are rotated into one another oppositely, for instance, if the {j, it}-plane 
is rotated by an angle a, and the {jt, i}-plane is rotated by the angle —a, then it turns out 
that the multiplication table is preserved. We have thus constructed a 1-parameter family of 
automorphisms. There are three ways of pairing up the three pairs of units in this way, but only 
two are independent. Since there are 7 different units that can be pointed to, the dimension of 
this group is again 14. 

In what follows, we will need not only G2, but also SU(3), the subgroup of G2 that fixes one 




Figure 4. A second class of elements of C?2- These transformations are not contained in the 
preferred SU (3) that fixes £. 



of the octonionic units. Since t is in the middle of our multiplication table, we will, without 
loss of generality, choose it to be the unit that is fixed. We see that the G2 transformation in 
Figure 3 fixes I and is therefore in SU(3), but a G2 transformation involving either of the other 
two pairs that point to k£ will not fix £. To be symmetric, we choose the linear combination 
of transformations shown in Figure 4 ({j, i£} and {j£, 1} both rotate by a and {k, £} rotates by 
—2a) to be the G2 transformation that points to k£ that is not in SU(3). If we choose to point in 
turn to each of the six units that are not £, we have six G2 transformations that are in SU(3) and 
six that are not. What about the remaining two G2 transformations? These are transformations 
that point to £. One such transformation, shown in Figure 5, rotates {i£, i} by a and {j£,j} by 
—a. There are three transformations of this type, all of which fix £ and are therefore elements 
of SU(3), but only two are linearly independent. Any two of these transformations complete 
the 8-dimensional Lie group SU(3). 

Yet another way to describe G2 is in terms of inner automorphisms, that is, transformations 
of the form 

x 1 — y axa^ 1 (21) 
Inner automorphisms always preserve an associative multiplication rule, since 

(axa~ 1 )(aya~ 1 ) = a(xy)a^ 1 (22) 



However, this condition is nontrivial over the octonions, since the parentheses cannot be moved. 
As shown in [9], (22) holds for all x,y G O if and only if a is a sixth root of unity. That is, the 
inner automorphisms of the octonions are given by (21) where 

a = e nwS/3 (23) 

where n G Z and s is any pure imaginary unit octonion. As further discussed in [9], any G2 
transformation can in fact be generated by a finite sequence of nested transformations of the 
form (21), with a given by (23). 



Figure 5. A third class of elements of Gi- These transformations are contained in the preferred 
517(3) that fixes I. 



5. The Jordan Eigenvalue Problem 

In previous work [5], we solved the Jordan eigenvalue problem, namely the eigenmatrix problem 

A o V = AV (24) 

where A and V are both 3x3 octonionic Hermitian matrices. Unlike the right eigenvalue problem 
Av = vX considered in [10], the Jordan eigenvalue problem (24) admits only real eigenvalues, 
which do solve the characteristic equation for A, namely [11] 

- det(A - AX) = X s - (tr A) X 2 + a(A) X - (det A) 1 = (25) 

where X denotes the identity matrix, o~(A) is defined by 

a(A) = i ((tr A) 2 - tr{A 2 )) = tx{A * A) (26) 
the operation * denotes the Freudenthal product 

A*B = AoB--(a tr(B) + B tr(X)) + ^ ( tr(^) tr(B) - tr(A o B)) X (27) 
and the determinant is defined unambiguously by 

det(A) = i tr ((A * A) o A) (28) 

The Jordan and Freudenthal products are generalizations of the standard dot and cross products. 

Just as in the more familiar complex case, normalized eigenmatrices for each nondegenerate 
eigenvalue are primitive idempotents, and the degenerate case can be handled using Gram- 
Schmidt orthogonalization. Furthermore, each primitive idempotent is in fact an element of the 
Cayley-Moufang plane OP 2 , which can be characterized as 

OP 2 := {V G H 3 (0) : V o V = V and tr V = 1} (29) 
= {V € H 3 (0) : V * V = and tr V = 1} 



It is straightforward to show from the first condition in (29) that the components of any element 
V G OP 2 must lie in some quaternionic subalgebra of O, which of course depends on V. Put 
differently, the associator of the (independent, off-diagonal) components of V, denoted [V], must 
vanish. But quaternionic primitive idempotents have (nonunique) "square roots", V = , so 
that we can also write 

OP 2 = : * G O 3 , [*] = 0, = 1} (30) 

where [\P] denotes the associator of the components of fy. We refer to such 3-component 
octonionic column vectors ^ as Cayley spinors. 

Putting this all together, any 3x3 octonionic Hermitian matrix A can be expressed as the 
sum of the squares of quaternionic columns, which are orthogonal under the Jordan product, 
that is 

3 

A = Y,Wi ( 31 ) 

1=1 

in terms of primitive idempotents Vj = G OP 2 and their corresponding eigenvalues Aj. 

6. The Structure of E 6 

The automorphism group of the Jordan product (2) (and consequently also of the Freudenthal 
product (27)) is the exceptional group F4, and the group which leaves the determinant (28) 
invariant is a particular real form of the exceptional group E$ . These groups can be interpreted 
as F4 = SU(3,0) and Eq = 5L(3,0), as we now show; for further details see [12]. 

In previous work [9], Manogue and Schray showed how to write the (double cover of the) 
Lorentz group 50(9, 1) as SL(2, O), with the action given by 

X .— > MXM 1 " (32) 

where X G H2(0), the 2x2 Hermitian matrices with octonionic components. The key to 
that work was to give an explicit set of basis transformations — the rotations and boosts in 
coordinate planes — which were compatible with the spinor representation in the sense that if 
9 G O 2 transforms like 

1 — ► M8 (33) 
then there are no associativity problems in the vector transformation 

M(6»6» t )M t = (M6»)(M6i) t (34) 

Any such basis transformation M G 5L(2,0) can be immediately reinterpreted as a 3 x 3 
transformation A4 via 

M 0' 



M = 1 (35) 



and it is straightforward to verify that any such A4 preserves the determinant of X G H%( 
and is therefore an element of E^- 

How many such transformations are there? We first give the basis transformations for the 
simpler case of 5L(2,C), adapted from [9] and rewritten as elements of Eq. When interpreting 
these transformations, it is helpful to recall that, in this case, 

X=( 4+ : r tV ) P6) 

We have the three rotations 

M xy = I c<°r2 j (: 5 7) 




M yz = 

and the three boosts 
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written in terms of the single imaginary unit i. To generate a 3 x 3 representation of 5L(2,0), 
start by replacing i in turn by each of the other imaginary units, yielding 18 new transformations, 
for a total of 24. The remaining 21 transformations in SL(2, O) are precisely the rotations of 
the imaginary units with each other, which make up an SO(7) (really Spin(7); we are being 
casual about double covers). As shown in [9], these rotations are obtained by nesting, that is by 
transformations of the form 

X i — > M.2 (MiXM^ M\ (41) 

where each corresponding M represents a "flip", that is, a pure imaginary multiple of the (2 x 2!) 
identity matrix. Thus, a typical M. takes the form 



M 




(42) 



where it is important to note that A4 is not a multiple of the (3 x 3) identity matrix. 

We are now ready to count the basis transformations of Eq. At first sight, it appears we 
have three copies of SL(2, O) — simply repeat the embedding (35) with the two other obvious 
block structures. We call these three copies type I, II, and III. However, this yields 3 x 45 = 135 
transformations, and, while these transformations do indeed generate all of Eq, it is clear that 
they can not be a basis, since the dimension of Eq is only 78. 

Let's try again. Each of these three copies of 5L(2, O) = SO(9, 1) contains a copy of SO(8). 
A famous property of SO(8) called triality asserts in this context that these three copies of 
SO(8) in fact consist of the same Eq transformations (but labeled differently), so we should 
count these copies only once. But 50(8) has 28 elements, to which we must add 3 copies of 
the 8 rotations needed to get to 5*0(9), then 3 copies of the 9 boosts needed to get to 50(9, 1), 
resulting in 28 + 3x8 + 3x9 = 79 transformations. A bit of thought reveals that the 3 copies 
of the tz-boost (39) are not independent; removing one of them correctly yields an explicit set 
of 78 basis transformations for Eq, also justifying the interpretation Eq = SL(3,0). 

It is worth pointing out that, due to triality, only the 14 G2 transformations need to be 
written in the nested form (41). Remarkably, the remaining 14 50(8) transformations can 
all be expressed using the type I transformation (37) and its types II and III variants. The 
former are just the usual 7 rotations needed to get from 50(7) to 50(8), but the latter yield 
an unnested description of the 7 non-G2 transformations in 50(7), which take the form 
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(43) 



Each of these transformations rotates 3 octonionic planes by the same amount, and can therefore 
be thought of "phase" transformation. 

What about F 4 ? Note that we have described 27 — 1 = 26 boosts, and 78 — 26 = 52 rotations. 
So our Eq is the real representation with 26 boosts, commonly written as £^(-26) , with the 
number in parentheses denoting the number of boosts minus the number of rotations. It is 
straightforward to show that F4 preserves the trace of elements of H^iO), corresponding to the 
timelike direction; this is the compact representation of F4, consisting precisely of the rotation 
subgroup of this real form of E%. These considerations justify the interpretation = SU (3, 0). 

Returning to the characteristic equation (25), not only does E$ preserve the determinant, it 
also preserves the condition a = 0. But these two coefficients control the number of nonzero 
eigenvalues — 3 if det.4 / 0, 2 if det.4 = / cr(A), and 1 if det.4 = = a (A). Thus, 
Eq preserves the number of nonzero eigenvalues of A, and hence the number of terms in the 
decomposition (31) with nonzero eigenvalue. 



7. Symmetry Breaking and Particle Physics 

In order to apply this formalism to elementary particle physics, we break the full symmetry and 
explore how a Jordan matrix transforms under various subgroups of Eq. One such symmetry 
breaking occurs when we choose a preferred SX(2,©) subgroup of SL(3, O), as in (35). This 
leads us to impose a block structure on H^(0), so that 

t) ^ 

where P € H2(0) transforms like a 10-dimensional momentum vector, ip £ O 2 transforms like a 
(Majorana-Weyl) spinor, and n € M is a scalar. Direct computation shows that 

V*V=[**l-f (45) 
\ (P^) f detPy/ 

where tilde denotes trace reversal, that is, P = P — tr(P) I. As shown in [14], the massless, 
momentum-space Dirac equation in 10 dimensions can be written 

PV> = (46) 

which implies the nonlinear constraint 



detP = (47) 

The general solution of (46) and (47) is of the form 

V> = 0£ (48) 

P = 99^ (49) 

where the components of 9 <G O 2 lie in the complex subalgebra of O determined by P and £ G O 
is arbitrary. Using (45), these equations are seen to be precisely the same as 

V * V = (50) 



where 



n=\i\ 2 



(51) 



Thus, (normalized) solutions of the Dirac equation are precisely elements of OP 2 , and therefore 
the squares of Cayley spinors. 

In previous work [14], we discussed solutions of the Dirac equation in the form (46). 
Remembering that solutions of (50) are quaternionic, and reducing 10 spacetime dimensions 
to 4 by the simple expedient of choosing a preferred complex subalgebra of O, we used spin 
eigenstates and particle/antiparticle projection operators as usual to identify particle states 
within our division algebra formalism. This procedure identified a spin-^ massive particle with 
two spin states, namely 




where the direction of the arrow indicates the z-component of the spin, and where the second 
equality in each case gives the momentum vector. Comparison with (36) shows that the x, y, 
and z components of the momentum vanish; these states are given at rest. Similarly, there is an 
analogous antiparticle with two spin states. The procedure also identified a left-handed massless 
particle, which when moving in the z-direction takes the form 

"* = (*) ^ t= (o !) ^ 

together with a massless particle of the opposite helicity, which when moving in the z-direction 
takes the form 

0,= (?) =(o ?) (55) 

Any other 2-component quaternionic column can be identified as an appropriately rotated and/or 
boosted superposition of these particles in the usual way. Noting that the first two particles 
carry an octonionic label (k), it is straightforward to generalize these particles to 3 generations 
of leptons, labeled by i, j, k; the remaining particle does not have an octonionic label. Thus, 
an octonionic description of the Dirac equation in 10 dimensions yields a particle spectrum 
containing precisely 3 generations of leptons, each with a single-helicity, massless neutrino, 
together with a single "sterile" neutrino of the opposite helicity but with no generation structure. 

8. Discussion 

We have shown how to break the symmetry group Eq so that a Lorentzian 3 + 1 dimensional 
momentum space emerges, together with internal symmetries that describe the correct 
spin/helicity transformations on Cayley spinors to describe leptons. Furthermore, precisely three 
generations of such leptons exist which respect the octonionic structure of the transformations. 
Contained naturally within this description of leptons and their symmetries are three massless 
left-handed neutrinos and a single, sterile, right-handed neutrino. 

We have learned several important lessons along the way. First, in spite of the non- 
commutativity and non-associativity of the octonions, everything that one might want to do can 
be made to work if one defines everything carefully. Second, when working with the octonions, 
it is important to make the Lie group structure primary, rather than the Lie algebra structure. 
Some of the nested group transformations described in (41) cannot be described in terms of the 
exponentials of any Lie algebra transformations. Finally, in order to keep a complex structure on 
the Lie algebra from interfering with the octonionic units in the matrices, it is important in the 
symmetry-breaking process to look at real Lie subalgebras rather than complexified ones. We 




A x = su('2.£) 

Figure 6. A map of Eq, taken from [13]. 

call the reader's attention to the recent work of Aaron Wangberg [13], in which the important 
real forms of the subgroups of E$ were identified, using a division-algebra perspective. A map 
of these subgroups, taken from [13], appears in Figure 6. 
We conclude with several speculative questions: 

(i) In this paper, we have described leptons in terms of "1-squares," i.e. the squares of Cayley 
spinors. But the most general octonionic hermitian matrix can be a linear combination of 
up to three such squares. Is it possible to describe the meson and baryon sectors of standard 
particle physics as "2-squares" and "3-squares," respectively? 

(ii) Does the group Eq, broken in the way we have discussed, describe both the standard model 
interactions and Lorentz transformations in 3 + 1-dimensions? 

(hi) After the symmetry breaking described in Section 7, some of the Lie group transformations 
of the form (21), (23) are no longer connected to the identity. Do these discrete 
transformations correspond to discrete conserved quantities such as charge? 
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